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Q\ • A major property of confined Fermi system is the quantisation of single particle 

motion. It leads to a bunching of levels in the single particle spectrum, known as 
shells, and gives rise to magic numbers in finite Fermi systems. Consequently, a 
spherical symmetry leads to very strong shell effects manifested in the stability of 
(*C) ■ the noble gases, nuclei and metallic clusters |], 0- Recently the sequence of magic 

numbers of a two-dimensional harmonic oscillator has been observed in the addition 
energy for a vertical quantum dot, i.e. in the energy needed to place the extra 
electron into the dot ||. In contrast to nuclei and metallic clusters the properties 
of quantum dots can be controlled by men. The main aim of this talk to discuss 
O ' the manifestation of shell effects in the mean field approach in different mesoscopic 
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systems, i.e. in nuclei, metallic clusters and quantum dots. 



C""-. . Periodic motion and higher multipoles in nuclei and metallic clusters. 

When a spherical shell is only partially filled, a breaking of spherical symmetry, 
resulting in an energy gain, can give rise to a deformed equilibrium shape. Super- 
and hyperdeformed nuclei are among the most fascinating examples where devia- 
tions from the spherical shape are a consequence of strong shell closures giving rise 
to largest level bunching (largest degeneracy or lowest level density). Nowadays it 
is recognized that the fine structure in the mass spectra between magic numbers in 
metallic clusters could be explained via symmetry breaking mechanisms similar to 



the situation in nuclear physics. It is therefore accepted that cluster deformations 
can exist, and it is actually confirmed at least for clusters with A < 40 either by the 
Clemenger-Nilsson (CN) model (introduced by Nilsson M in nuclear physics and 
applied by Clemenger || for clusters) or by a self-consistent Kohn-Sham density- 
functional method (KS) with deformed jellium backgrounds [[7|~@- The need for 
multipole deformations higher than the quadrupole in the mean field approach has 
been recognized in nuclei and in metallic clusters in numerous calculations to ex- 
plain experimental data. For instance, the octupole deformed shapes still constitute 
an intriguing problem of nuclear structure, experimental as well as theoretical (see 
for review The hexadecapole deformation is essential for the understanding of 
equilibrium shapes and the fission process of super- and hyperdeformed nuclei . 



In the case of metallic clusters, the axial hexadecapole deformation is important for 



the interpretation of experimental data in simple metals |pj]| . Inclusion of higher 
multipoles leads, however, to a nonitegrable problem. In fact, the single particle 
motion turns out to be chaotic. It is expected that an increasing strength of the 
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corresponding multipole deformation increases the amount of chaos in the classical 
potential as quantified by an increasing Lyapunov exponent. Accordingly, in the 
transition ^from ordered to chaotic motion the quantum numbers lose their signifi- 



cance, and the system behaves like a viscous fluid [12]. Therefore, a disappearance 
of shell structure should be expected in the analogous quantum case. However, re- 
cently the occurrence of shell structure has been reported for many body systems 
like nuclei and metallic clusters O, inj at strong octupole deformation. A major 
conclusion of ]14j is that, albeit nonintegrable, an octupole admixture to quadrupole 
oscillator potentials leads, for some values of the octupole strength, to a shell struc- 
ture similar to a plain but more deformed quadrupole potential. This result shows 
that there is a tendency of the system to restore the original symmetry (degeneracy 
of levels) which is destroyed by the octupole term. 

The equilibrium deformation is ultimately related to the behavior of the single 
particle level density of the quantum Hamiltonian. According to the semiclassical 
theory |15| the frequencies in the level density oscillations of single particle spectra 
are determined by the corresponding periods of classical closed orbits. The short 
periodic orbits determine the gross shell structure, whereas contribution of longer 
orbits give finer details. Analysis of shell structure phenomena in nuclei in terms of 
classical orbits was started by the work of Balian and Bloch []I6] who studied the 
density of eigenmodes in a spherical cavity with reflecting walls. The importance of 
deformed shapes for nuclei led to the generalization of this analysis in considering 
a deformed ellipsoidal well and a deformed harmonic oscillator [|], |I7| . Since these 
problems are integrable, attention was given to regular motion. One of the first 
attempts to analyze shell structure phenomena and chaotic motion in a quadrupole 
deformed diffuse cavity was given in |Tj|. Our approach is based on the connection 
between shell structure phenomena in the quantum spectrum and ordered motion 
in the classical analogous case. As it was mentioned above, shell structures in the 
quantum mechanical spectrum are related to periodic orbits of the corresponding 



classical problem [!]],[ 16], [17| . The periodic orbits are associated with invariant tori 



of the Poincare sections. If the classical problem is chaotic, the invariant tori disinte- 
grate or disappear |25[] , and the shell structure of the quantum spectrum is affected 
by the degree of chaos [rS|-[23j|. The classical treatment is based on the secular 



perturbation theory |2£| and is particularly effective for a two degrees of freedom 
system. The technique uses action-angle variables of the unperturbed Hamiltonian 
and averages over the faster phase. Usually, prior to such an operation, a canon- 
ical transformation is necessary in order to remove the initial resonance from the 
unperturbed Hamiltonian (in our case the axial harmonic oscillator without higher 
multipoles). In the new rotating frame, one of the phases will only measure the 
slow variation of the variables about the original resonance which now becomes a 
fixed elliptic point. The problem is then treated by averaging over the remaining 
faster phase. For super- or hyperdeformed system the frequency of the oscillation 
in the ^-coordinate is twice as fast or faster than that in the z-coordinate. Since 
the Hamilton function is periodic in both angles, this approximation amounts to 
keeping the zero order term of its Fourier expansion in the fast moving angle. The 
integration of the Hamilton function over the angles can be done analytically, which 
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makes the approximation particularly attractive. The full Hamilton function, being 
a nonintegrable problem, is approximated by an effective Hamilton function that is 
obtained by averaging over the fast oscillating angle $ Q for the case considered. 

An axial harmonic oscillator Hamiltonian (AHO) which is deformed by arbitrary 
multipoles has the form 

H = ±(fi+£ + *)+U{r,*) (1) 

where 

U(r,#) = ^f- (V + yi + r 2 (A 3 P 3 (costf) + A 4 P 4 (costf) + ...)) (2) 

with r 2 = g 2 + z 2 , cos$ = z/r. The z-component of the angular momentum is 
denoted by p<p and Pk{cos$) is the k— th order Legendre polynomial. The parameter 
b characterizes an oblate and a prolate shape for < b < 1 and b > 1, respectively. 
For a non- vanishing A3, A4, . . . the problem becomes a nonintegrable two degrees of 
freedom system. Note that the only constant of motion is the total energy. Our 
interest is focused on the contribution of the octupole and hexadecapole deformation 
in super- and hyperdeformed system. Therefore, the general potential Eq.([|) is 
reduced to the form 

m 2 2 z 2 2z 3 - 3zg 2 8z 4 - 2Az 2 g 2 + 3g\ 

v ^ z) = ^( Q+ b 2+ X '^FTW + z^T? } - (3) 

The terms multiplied by A 3 and A 4 give rise to octupole and hexadecapole deforma- 
tions, the respective terms are proportional to r 2 P 3 (cos $) and r 2 P 4 (cos 

The axial symmetry of the potential given in Eq.(^) guarantees conservation of 
the ^-component of the angular momentum denoted by p^, the discussion will be 
focussed on p^ = unless indicated otherwise. After averaging and rewriting the 
action variables and the remaining angle d z in terms of the original momentum and 
coordinate values we obtain [^] 

av = 2m IT 9 cS ^ ^ 



where 



= + ^P^{4mA) - 3^F,(i 1 2; 4)) 
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A 4 (-£ 2 - 27z 2 + ' )]. (5) 

2 8 T.&T^ n 



and £ 2 = 2J e /(mu) = 2E Q / {muj 2 ) which is a constant within the approximation. 
Here K is the first elliptic integral. Note that the approximated Hamilton function 
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is separable in the two coordinates with the ^-motion being unperturbed. Conse- 
quently the frequencies are u e = to and u z = 2n/T z with 

, /"^max H V 

T Z = V2^ (6) 
J*** ^E-E e -U cS {z) 

The motion in the ^-coordinate is different from the unperturbed motion and its 
frequency depends on £ e , i.e. on the amount of the energy residing in the ^-motion. 
We obtain for the winding number 

^ = I( I + I ) (7) 

Uz 2 ^Ji/ h 2 Z 2A 3 + 8A 4 + 2A 3 + 8A 4 

This result is exact for H av when £ e = (all energy resides in the z-motion), but by 
the previous argument it applies for large part of phase space. Moreover, it serves 
as a useful and reliable guideline for the full problem. 

The oblate case {b < 1) does not lend itself to the same approximation procedure. 
As was reported in an octupole addition to an oblate quadrupole potential gives 
no contribution to the zeroth order when averaging over d z . If only a hexadecapole 
term is added, the same procedure can be applied and we obtain 

W eS (g) = + A 4 (4^L - 32^ 2 + 4^ 2 )] (8) 
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Q 2 + g 



with (1 = 2bE z /(muj 2 ). In this case, we find for £ z = 0, for the winding number 

CO 



-^ = bJl + 3X 4 . (9) 

uj z v 

Note that the admixture of the hexadecapole term to the oblate potential yields, 
quantum mechanically, an effective plain oblate case but with less deformation. 
This is in contrast to adding an octupole to a prolate potential where the effective 
deformation is increased. 

The quantum-mechanical results are in line with the classical predictions. We 
calculated the energy levels of the Hamiltonian 

H = H + X 3 H 3 + A 4 # 4 (10) 

where the diagonal matrix Ho comprises the axial harmonic oscillator, and H3 and 
i?4 are the octupole and hexadecapole terms, respectively. The consistent calculation 
of the matrix elements in a truncated basis is described in H . 

As a quantitative measure for shell structure we use the Strutinsky type analysis 



(see for details [H], Q). In accordance with the discussion above, the oblate su- 
perdeformed potential produces chaos when the octupole term is switched on. The 
quantum spectrum has the level statistics ascribed to quantum chaos [ 2D| . The hex- 



adecapole term alone does gives rise to new shell structure; for instance, for b = 2/5 
and A4 ~ 0.1 oblate super deformation, i.e. u) q /uj z = 1/2 was established. This is 



close to the value given by Eq.(0), the difference is explained in [22 . 
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There is remarkable agreement between the manifestation of shell structure for 
values of the strength parameters A3, A4, which coincide with the ones predicted by 
the classical perturbative approach and give rise to stability islands on the Poincare 
surfaces of sections. According to the classical approach only even multipoles can 
decouple the potential of Eq.(l) for an oblate deformation. Therefore, shell struc- 
ture can be supported for strongly oblate deformed nuclei or clusters only by even 
multipoles. A further interesting result is the mutual cancellation of the octupole 
and hexadecapole contribution in the quadrupole deformed system. The classical 
approach allows to determine the range of parameters of A3 and A4 where the cor- 
responding quantum mechanical problem of the quadrupole + octupole + hexade- 
capole potential yields shell structure resembling the plain quadrupole deformation 
p4|j . Along this curve the octupole deformation tends to increase the effective pro- 
late deformation whereas the hexadecapole term produces the opposite effect. In 
fact, the hexadecapole term can stabilize the octupole deformation in superdeformed 
systems, since the cancellation curve attains values of A3 which are larger than its 
critical value for A4 = 0. We may speculate that prolate superdeformed nuclei with 
rather strong octupole deformation could therefore exist. Finally we comment that 
even though the quantum mechanical treatment shows a certain degree of suppres- 
sion of classical chaos, the occurrence of a new shell structure which differs from the 
unperturbed case is clearly brought about by the nonlinear character of the problem. 

Filled shells in quantum dots. We choose the harmonic oscillator potential 
as the effective mean field for the electrons in an isolated quantum dot. While the 
electron-electron interaction is important for the explanation of certain ground state 
properties like special values of angular momenta of a quantum dot in a magnetic 
field |S7| , a smooth and finite potential which admits bound states can for the lowest 



few levels always be approximated by the harmonic oscillator potential Pq| . This 
has been confirmed by direct determination of a new effective oscillator frequency 
for two interacting electrons in an external parabolic potential |29| and by calcula- 



tions of the effective single particle levels within the density-functional theory for 
electron numbers N ~ 100 For small dot size and small number of electrons the 
confinement energy becomes prevalent over the Coulomb energy. The experimen- 
tal identification of the magic numbers of the two-dimensional harmonic oscillator 
without a magnetic field || is the most convincing argument for the validity of this 
approximation. The effect of an external homogeneous magnetic field can be calcu- 
lated exactly for a three dimensional (3D) harmonic oscillator potential irrespective 



of the direction of the field ||31|| . Our discussion here is based upon the 2D version 
of the Hamiltonian |31| including spin degree of freedom. The magnetic field acts 
perpendicular to the plane of motion, i.e. H = J2f=i hj with 
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o ^ (p- -Ay + — (^V + ^V) + ^a z B. (11) 

where A = [fxB]/2, B = (0, 0, B) and a z is the Pauli matrix. We do not take into 
account the effect of finite temperature; this is appropriate for experiments which 
are performed at temperatures kT <C A with A being the mean level spacing. In 
the following we use meV for the energy and Tesla for the magnetic field strength. 
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The effective mass which determines the orbital magnetic moment for the electrons 
is chosen as m* = 0.067m e . It leads to /iff = 15/is while the effective spin magnetic 
moment is fi* = 0.5/Xb. 

Shell structure occurs whenever the ratio of the two eigenmodes Q± of the Hamil- 
tonian flTT| ) (see Ref.|JT|) is a rational number with a small numerator and denom- 
inator. If we start with a circular dot (u x = LU y ), the shell structure is particularly 
pronounced if the ratio is equal to one (for B = 0) or two (for B pa 1.23 ) or three 
(for B pa 2.01 ) and lesser pronounced if the ratio is 3/2 (for B = 0.72) or 5/2 (for 
B = 1.65). The values given here for B depend on m* and w x>y . As a consequence, 
a material with an even smaller effective mass m* would show these effects for a 
correspondingly smaller magnetic field. The magic numbers (including spin) turn 
out, for B — 0, to be the usual sequence of the two dimensional isotropic oscilla- 
tor, that is 2, 6, 12, 20, . . . ||. For B pa 1.23 we find a new shell structure as z/the 
confining potential would be a deformed harmonic oscillator without magnetic field. 
The magic numbers are 2, 4, 8, 12, 18, 24, . . . which are just the numbers obtained 
from the two dimensional oscillator with u > = 2uj < (o;> and uj < denote the larger 
and smaller value of the two frequencies). Similarly, we get for B pa 2.01 the magic 
numbers 2, 4, 6, 10, 14, 18, 24, . . . which corresponds to u > = 3u> < . If we start from 
the outset with a deformed mean field, i.e. if we choose, say, uj x — (1 — f3)ui y with 
(3 > two major effects are found: (i) the degeneracies (shell structure) are lifted 
at B — depending on the actual value of (3, and (ii) the values for B at which 
the new shell structures occur are shifted to lower values. The significance of this 
finding lies in the restoration of shell structures by the magnetic field in an isolated 
quantum dot that does not give rise to magic numbers at zero field strength due to 
deformation. We mention that the choice (3 = 0.5 would shift the pattern found at 
B pa 1.23 to the value B = 0. It is the shell structure caused by the effective mean 
field which produces the maxima that are observed experimentally in the addition 
energy fi(A + 1) — fi(A) = Ea+i — Ea + e 2 /C 0. Here Ea is the single particle 
energy of the effective mean field in quantum dots, e 2 /C is the electrostatic energy 
and fi(A) is the chemical potential. The electrostatic energy is much larger than the 
difference Ea+i — Ea, however, it is the fluctuations (shell effects) of the difference 
that matters, at least for small quantum dots. The effect is similar to shell phenom- 
ena in nuclear physics and for metallic clusters. The analogy goes further in that, in 
an isolated small quantum dot, the external magnetic field acts like the rotation on 
a nucleus thus creating new shell structure; in this way superdeformation (axis ratio 
2:1) has been established for rotating nuclei owing to the shell gaps in the single 
particle spectrum. We now focus on the special cases which give rise to pronounced 
shell structure, that is when the ratio Q + /Q- = k = 1, 2, 3, . . ., and analyse in detail 
the circular shape (u> x = u y = We find for the magnetization 

M = /if (1 - ^)(£_ " k T,+) ~^<S Z > (12) 

with E± = Ef(n± + 1/2) j 0, lj l = aB = ^-B and = ha being the effective 
Bohr magneton using the effective mass m*. The eigenmodes are Q± — (Q ± ojl) 
with Q = Ju) 2 + io\ [[32] . 
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For completely filled shell < S z >= 0, since, for the magnetic field strenghts 
considered here, the spin orientations cancel each other. ^From the orbital motion 
we obtain for the susceptibility 

eff 2 

^/^-MNE^EJ da) 

It follows from Eq.(|H|) that, for a completely filled shell, the magnetization owing to 
the orbital motion leads to diamagnetic behaviour. For zero magnetic field (k = 1) 
the system is paramagnetic and the magnetization vanishes (J2- — The value 

k = 2 is attained at B m 1.23. When calculating J2- an d J2+ we now have to 
distinguish between the cases, where the shell number N of the last filled shell is 
even or odd. For the last filled shell number even one finds 

E + = ^(iV + 2)[(iV + 2) 2 + 2] (14) 

J2_ = ^(N + l)(N + 2)(N + 3) (15) 

and M = -fj,f{L-u L /Q)(N + 2)/2. 

For the last filled shell number odd one finds 

E + = = + + 2)(N + 3) (16) 

and M — 0. Therefore, if f2 + /fi_ = 2, the orbital magnetization vanishes for the 
magic numbers 4, 12, 24, . . . while it leads to diamagnetism for the magic numbers 
2, 8, 18, .... A similar picture is obtained for Q + /Q— = 3 which happens at B « 2.01: 
for each third filled shell number (magic numbers 6, 18, . . .) the magnetization is 
zero. Since the results presented are due to shell effects, they do not depend on 
the assumption u x /uj y = 1 which was made to facilitate the discussion. The crucial 
point is the relation Q + /Q- = k = 1,2,3, .. . which can be obtained for a variety 
of combinations of the magnetic field strength and the ratio uj x /uj y . Whenever the 
appropriate combination of field strength and deformation is chosen to yield, say, 
k = 2, our findings apply. 

In conclusion: our analysis for nuclei, metallic clusters and quantum dots is 
based on a rather simple model, and it is true that a realistic single particle spec- 
trum is poorly reproduced by the the harmonic oscillator. However, the lucidity and 
transparence of the model in describing the phenomena of extreme deformations for 
nuclei and metallic clusters is superior to any realistic model. The shell structure 
effects observed for the addition energy of a small isolated quantum dot provide 
a reliable basis for use the 2D version of the harmonic oscillator Hamiltonian. At 
certain values of the magnetic field strength shell structures appear in a quantum 
dot, also in cases where deformation does not give rise to magic numbers at zero 
field strength. Measurements of the magnetic susceptibility are expected to reflect 
the properties of the single particle spectrum and should display characteristic pat- 
terns depending on the particle number. At certain values of the magnetic field 
and electron numbers the orbital magnetization vanishes due to shell closure in the 
quantum dot. This property could be of interest in applications because it enables 
control of the electron number in small isolated quantum dots. 
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